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It is shown that the copropagating three-wave-mixing parametric process, with ap- 
propriate type-II extended phase matching and pumped with a short second-harmonic 
pulse, can perform spectral phase conjugation and parametric amplification, which 
shows a threshold behavior analogous to backward wave oscillation. The process is also 
analyzed in the Heiscnbcrg picture, which predicts a spontaneous parametric down 
conversion rate in agreement with the experimental result reported by Kuzucu et al. [Phys. 
Rev. Lett. 94, 083601 (2005)]. Applications in optical communications, signal processing, 
and quantum information processing can be envisaged. © 2008 Optical Society of America 
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1. Introduction 

In contrast with the more conventional optical phase conjugation schemes that perform phase conjugation 
with spectral inversion, 1 spectral phase conjugation (SPC) is the phase conjugation of an optical signal in the 
frequency domain without spectral inversion. Equivalently, in the time domain, SPC is the phase conjugation 
and time reversal of the signal complex pulse envelope. 2 SPC is useful for all-order dispersion and nonlinear- 
ity compensation, 3,4 as well as optical signal processing. 5 Although SPC has been experimentally demon- 
strated using photon echo, 6 ' 7 spectral hole burning, 8,9 temporal holography, 3 spectral holography, 10 and 
spectral three-wave mixing (TWM), 11 all the demonstrated schemes suffer from the use of cryogenic setups, 
non-realtime operation, or extremely high pump energy. Pulsed TWM 12 and four-wave-mixing (FWM) 2, 13 
processes in the transverse-pumping geometry have been theoretically proposed to efficiently perform SPC, 
but have not yet been experimentally realized. All the holographic and wave-mixing schemes also have strict 
requirements on the transverse beam profile of the signal, limiting their appeal for simultaneous diffraction 
and dispersion compensation applications. 

There is a correspondence between classical SPC and quantum coincident frequency entanglement, as 
shown in Ref. 14 for the transversely pumped TWM 12, 15 and FWM 2, 13 processes. It is then interesting to 
see if other coincident frequency entanglement schemes arc also capable of performing SPC, when an input 
signal is present. This paper studies one of such schemes, which makes use of extended phase matching 
(EPM) 16 and has been experimentally demonstrated 17 in a periodically-poled potassium titanyl phosphate 
(PPKTP) crystal. 18 It is shown in Section [3] for the first time to the author's knowledge, that this EPM 
scheme is indeed capable of performing SPC and optical parametric amplification (OPA), more efficiently 
than previous proposals. 

The analysis also yields a surprising result, namely that the parametric gain can be theoretically infinite 
even for a pump pulse with finite energy, analogous to backward wave oscillation, where countcrpropagating 
waves are parametrically coupled and can give rise to mirrorless optical parametric oscillation (OPO). 19-25 
The reason for the similarity is that, in the scheme presented here, even though the signal and the idler 
copropagate with the pump pulse in the laboratory frame, they counterpropagate in the frame of the moving 
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pump pulse, because one is faster than the pump and one is slower. Hence the moving pump pulse provides 
both an effective cavity and parametric gain, leading to oscillation. In reality, however, the interaction 
among the pulses should be ultimately limited by the finite device length. It is shown in Section 0] with 
a Laplace analysis, that the parametric gain should abruptly increase above the threshold, where infinite 
gain is predicted by the Fourier analysis, but a finite medium length would always limit the gain to a finite 
value. Still, as previous proposals of TWM mirrorlcss OPO have never been experimentally achieved due to 
the requirement of a continuous-wave (CW) pump and the difficulty in phase matching counterpropagating 
waves, the presented analysis suggests the exciting possibility that mirrorless OPO can be realized with an 
ultrashort pump pulse and a practical poling period for phase matching of copropagating modes, if a long 
enough medium can be fabricated and parasitic effects can be controlled. By analyzing the scheme in the 
Heisenberg picture in Section [5] a high spontaneous parametric down conversion rate is also predicted, in 
excellent agreement with the experimental result reported in Ref. 17. The result should be useful for many 
quantum information processing applications, such as quantum-enhanced synchronization 26 and multiphoton 
entanglement for quantum cryptography. 27 Finally, numerical results are presented in Section which 
confirm the theoretical predictions. 

2. Setup 




Fig. 1. Schematic of spectral phase conjugation (SPC) via type-II extended phase matching 
(EPM). The signal and idler pulses, in orthogonal polarizations, have carrier frequencies 
of uj s and u>i, while the pump pulse has a carrier frequency of lo p = lu s + uii. The EPM 
condition requires that the signal and the idler counterpropagate with respect to the pump, 
which should be much shorter than the input signal. 



Consider the copropagating TWM process (Fig.QJ, assuming that the basic type-II phase matching condi- 
tion (k s + fcj = k p + 27r/A), with a quasi-phase- matching period A, is satisfied. The coupled-mode equations 
are 

dA p , , dA 



dz ' k ' p ~dt~ ~ ^"P^-sAi, (1) 
^ + K ^ t=jx , Ar A h (2) 

?£ + k f-«L = - iXlA;A „ (3 ) 

where A p is the pump pulse envelope of carrier frequency lo p , A s j are the signal and idler envelopes of 
frequency uj s and cui respectively, k' psi are the group delays of the three modes, Xp,s,i = i^p,s,iX^ / (^ cn p,s,i) 
are the nonlinear coupling coefficients, uj PiS ,i are the center frequencies of the modes such that lo s + w% = u> p , 
and n PiSj i are the refractive indices. Group- velocity dispersion within each mode and diffraction are neglected. 
Define r = t — k'z as the retarded time coordinate that follows the propagating pump pulse. The change of 
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coordinates yields 



dA p 



dz 
dA* 



dr 
dA* 



jX P A s Ai, 
JXsApA* , 
-j Xi A*A s . 



(4) 
(5) 
(G) 



dz v 1 p/ dr 

Throughout the theoretical analysis, the pump is assumed to be undepleted and unchirped, so that A p = 
A p o(t — k'z) = A p o(t), hereafter regarded as real without loss of generality. 

3. Fourier Analysis 

Equations JSJ and © are space- invariant, if the nonlinear medium length L is much longer than the signal 
or idler spatial pulse width in the frame of z and r, or 
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T, 
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where T s .i is the signal or idler pulse width. One can then perform Fourier transform on the equations with 
respect to z, as defined by the following, 



A s {k,t) = 
A*(k,t) = 



A s (z, t) exp(— jnz)dz, 
A*(z, t) exp(— jnz)dz. 



(8) 
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Notice that A* is defined as the Fourier transform after the conjugation of A4. The coupled-mode equations 
become 



dA 

j K A s + (k' s - k' p )-g^ = 3XsA p o{t)A* 
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Consider the case in which j s and 7, are non-zero and have opposite signs, implying that the signal and 
the idler propagate in opposite directions with respect to the pump. This can be achieved for a range of 
wavelengths in KTP. Without loss of generality, assume that 7 S > and < 0, so that k' s > k' p > k[. 
Making the following substitutions, 



A = ^/rA s exp(j— r), B = A* exp(j— t), 

Is li 



one obtains 



dA 
dB 



= J 
= 3 
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Due to linear space invariance, the wave-mixing process cannot generate new spatial frequencies (k) for A 
and B. The magnitude of k then depends only on the initial bandwidths of A and B, and is on the order 
of 2Tt r y s _i/T s j. As a result, if the pump pulse width T p is much shorter than the minimum period of the 
detuning factor exp[±j«(l/7 s — 1/7j)t], or 



T p « 



2tt 



k(1/7s - I Hi) 7a,i(V7s - !/7i) 



(16) 



the pump can effectively sample the detuning factor, say, at r = 0. Defining a normalized coupling function 

9(r) = 



XsXi 



Isli 



^po(t), 



two simple coupled-mode equations are obtained, 



— =jg( T )A. 
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Because the signal and the idler countcrpropagate with respect to the pump, the signal should begin to mix 
with the pump at the leading edge of the pump pulse, say at r = — T p /2, while the idler should begin to mix 
at the trailing edge of the pump, say at r = T p /2. The solutions of Eqs. (|18|) and (|19|) can then be written 
as 



A(k,t) = sec(G)< A(k, 
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The input signal pulse is required to be placed in advance of the pump (by t s >> T s ), and the input idler 
pulse to be placed behind the pump (delayed by ij >> Tj), so that the signal and the idler only overlap the 
pump pulse inside the nonlinear medium. Consequently, the output solutions are 



A s (L,t) = A s0 {t - k' s L + t s ) sec(G) + 7-7=^, ( - -(* - K L - tan(G), 

y/r r 

Ai(L,t) = Aio{t-kiL-ti)sec(G) + jVrA* (- r(t - k\L + t s )) tan(G). 



(23) 
(24) 



To see how the device is able to perform SPC, assume that the center frequencies of the two modes are the 
same, u> s = u>i, Xs — Xi-, ancl the typc-II EPM condition, 



k s + k i — 2k p , k s 7^ fc i7 



(25) 



which depends on the material dispersion properties and typically occurs at a single set of center frequencies, 
is satisfied. 16 Then r = 1, and the output idler becomes the phase-conjugated and time- reversed replica of 
the input signal, if the input idler is zero. SPC is hence performed. The SPC efficiency 77, or the idler gain, 
defined as the output idler fluence divided by the input signal fluence, is 



C x \M0,t)\Mt 



tan 2 (G). 



(26) 
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This SPC efficiency can be fundamentally higher than that of the transversely pumped TWM device 12 due 
to two reasons. One is the copropagation of the three pulses, which makes G higher than a similar parameter 
in the latter case by a factor of (1 — k^/k^)" 1 , on the order of 40 for KTP. The second reason is that for 
rj > 1, due to the tangent function dependence, the SPC efficiency of the EPM scheme increases with respect 
to G much faster than that of the latter, which only depends on a similar parameter exponentially That 
said, the transversely pumped FWM device 13 can still be more efficient in the small gain regime t] < 1 if a 
highly nonlinear material, such as polydiacetylene, is used. Furthermore, the EPM device requires a longer 
nonlinear medium length by a factor of (1 — k'/k' s ) , and depends crucially on the material dispersion, thus 
severely limiting the flexibility in the choice of operating wavelengths. 

Equations (|23|l and (|24|l are obtained from the analysis of the coupled- mode equations J5j and JBJ), after 
Fourier transform with respect to z is performed. The solutions are therefore formally valid only when 
the nonlinear medium length L goes to infinity. In practice, in the moderate gain regime rj ~ O(l), the 
approximation given by Eq. (JZJ) should be adequate, where the length L can be, say, ten times larger than 
the signal spatial pulse width in the frame of z and r. Numerical analysis in Section [S] will validate the 
accuracy of the Fourier solutions. 

4. Laplace Analysis 

Intriguingly, the Fourier solutions, Eqs. i|23|) and Q24[l. have the same form as those of backward wave 
oscillation, 19-25 suggesting that the device studied here, with an ultrashort pump pulse and a practical 
quasi-phase-matching period (A = 46 /im as reported in Ref. 18), can also perform mirrorlcss OPO, as long 
as k' s i > k' p > k[ s . However, the prediction of infinite gain is based on the assumption of infinite medium 
length and therefore may not be valid. In this case, Laplace transform should be used. 

For the CW-pumped mirrorless OPO schemes, a Laplace analysis 28 with respect to time shows that beyond 
threshold, poles appear on the right-hand plane in the Laplace domain, meaning that the temporal impulse 
response increases exponentially with time, leading to self-oscillation when enough time is elapsed. The 
same procedures of utilizing the two-sided Laplace transform 29 as in Ref. 28 are followed here in order to 
be consistent with the relevant literature, but since the proposed scheme is the opposite limit of the CW 
devices, the Laplace transform should be performed with respect to z instead, 



A s (p,t)= A s (z,T)exp(-pz)dz, (27) 

J —OO 

A*Cp,r)= / A*(z,T)exp(-pz)dz. (28) 



For simplicity but without affecting the qualitative behavior of the solutions, it is assumed that the pump 
pulse is square, there is no input idler, 7 = 7 S = — 7i, and x = Xs = Xi- The output solutions in the Laplace 
domain are then given by 



A ^~Y > = P + VT^ZiGVT^) 1 ^-^ <»> 

p =—J—> G = xApo( — )- (31) 
xApo 7 

If we let p = jn, the transfer functions in Eqs. l|29|l and 1|30[) are well-known to be low-pass filters, 30 the 
bandwidth of which decreases as G increases. If the spatial bandwidth of the input signal, on the order of 
7/T s , is much smaller than the bandwidth of the low-pass filters, the transfer functions can be regarded 
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as flat-top functions, and by plugging P = in Eqs. (|29[l and (|3L)[) . the Fourier solutions in Eqs. 12-il) and 
(12411 are recovered. For G « 1, the transfer functions are sine functions with a bandwidth ~ l/T p , so the 
Fourier solutions are valid if T p << T Sl which is essentially the same assumption used in the Fourier analysis, 
Eq. As G increases and the filter bandwidth decreases, however, the Fourier solutions become less and 
less accurate for a finite-bandwidth input signal. 

The poles of the transfer functions, p^, can be obtained by setting the denominator of Eqs. (|29fl and (|30|) 
to zero, 

= 0. (32) 



Poo + J(xApo) 2 -P^cot GJl-p^/ixApo) 



Positive Poles of Transfer Functions 
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Fig. 2. Normalized poles Poo/{xA p o) plotted against G, obtained by numerically solving 
Eq. (|32|l . indicating the onset of spatial instability beyond the threshold G > ir/2. More 
poles appear as G is increased. 

Figure |2l plots the normalized poles Poo/(xA p o) against G. Positive poles begin to appear when G > ir/2, 
hence the spatial impulse response increases exponentially with respect to z beyond threshold. 

It is interesting to compare the scheme studied here with the case in which the pump, signal and idler 
have degenerate group delays (k' p = k' s = k^). 31 The coupled-mode equations of the latter case are 

~~ f g~ ~ = jxApo(T)Al s (z, t), (33) 

where the r derivatives vanish. The solutions are easily seen to be 

A s , t (z,T) = A Sil (0,T)cosh[xA p0 (T)z} + jA* s (0, t) sinh[xA p0 (r)z]. (34) 

This corresponds to the G — ► oo limit of the former scheme, where Poo/ (xA p o) — ► 1 and all the poles approach 
the growth rate of the degenerate case, xA P o- 

5. Spontaneous Parametric Down Conversion 

Given the input-output signal- idler relationship in Eqs. I|23|l and (j24(l . it is straightforward to obtain a 
quantum picture of the parametric process in the moderate gain regime by replacing the signal and idler 
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envelopes with Heisenberg operators, so that 



A s = A s0 sec(G) + jA\ Q tan(G), 
Ai = jA f s0 tan(G) + A M sec(G). 



(35) 
(36) 



If the inputs are Fock states, 



n Sii ee (Al.A s>i ) = (A^Alj - 1, 
(Al A i0 ) = (A\ A s0 ) = (A s0 Al ) = (A i0 Al ) = o 



(37) 
(38) 



The average output photon number of each mode is 



n s = n s0 sec 2 (G) + (n i0 + 1) tan 2 (G), 
rii = n i0 scc 2 (G) + {n s0 + 1) tan 2 (G). 



(39) 
(40) 



The average number of spontaneously generated photon pairs per pump pulse is therefore the same as the 
idler gain, or i] = tan 2 (G). Moreover, the unitary transform given by Eqs. i|35|) and H36|) has the same form 
as the CW FWM process. One then expects the photon wavcfunction to be similarly given by 32 



where \n) S: i is the Fock state in the signal or idler mode. The scheme thus has a significant advantage in 
efficiency and robustness for multiphoton entanglement, compared with other schemes that often require 
feedback. 33 The efficient multiphoton coincident frequency entanglement should be useful for quantum- 
enhanced synchronization 26 and quantum cryptography applications. 27 

The preceding quantum analysis assumes that there is only one spatial mode in each signal or idler mode, 
and is accurate only when the Fourier solutions are accurate. This restricts the applicability of the quantum 
analysis to the moderate gain regime 77 ~ O(l), depending on how closely the assumption in Eq. (0 is 
observed. It is beyond the scope of this paper to investigate what happens in the quantum picture when 
more than one spatial modes are involved, but qualitatively, one expects that each spatial mode should 
have a varying parametric gain depending on the spatial frequency, as suggested by the Laplace solutions in 
Eqs. 1|29[) and l|3U|l . so the photon wavefunction would be given by a superposition of simultaneous eigenstates 
of spatial frequency and photon number. 

Using the parameters described in Refs. 17 and 18, where Ao = 1584 ran, y^ 1 = 7.3 pm/V, no = 2, 
7 = 1.5 x 1CT 10 s/m, T. p = 100 fs, average pump power = 350 mW, diameter = 200 ^m, and pump 
repetition rate f rep = 80 MHz, the spontaneously generated photon pairs per second is theoretically given 
by frep tan 2 (G) f=a frepG 2 = 3.6 x 10 6 /s, in excellent agreement with the experimental result reported in 
Ref. 17, which is ~ 4 x 10 6 /s. G is then given by ~ 0.2, so the operations of SPC, OPA, and multiphoton 
entanglement (G > tt/4) should be realizable by increasing the pump field amplitude. 

6. Numerical Analysis 

Equations JSJ and © are solved numerically via a Fourier split-step approach to confirm the above theoretical 
predictions. Fig. [3] plots the intensities and phases of the input signal, output signal, and output idler from 
the numerical analysis when G = 7r/4. The plots clearly show that the output idler is the time- reversed and 
phase-conjugated replica of the signal. 

Figure 01 plots the numerical signal gain and idler gain compared with Fourier theory for < G < 7r/3. 
The numerical results arc all within 3% of the theoretical values. 



00 




(41) 



n=0 



7 




Fig. 3. Plots of intensity and phase of input signal, output signal and output idler, from 
numerical analysis of Eqs. JHJ) and ©. Parameters used are k' p = 1/(1.5 x lC^ms^ 1 ), k' s = 
1.025A£, ki = 0.975/fc^, T p = 100 fs, T s = 2 ps, L = 10 cm, t s = 4T S , beam diameter 
= 200 Mm, A s0 = 0.5exph(t - 2T S ) 2 /(2T S 2 )] - exp[-(l + 0.5j)(t + 2T S ) 2 /(2T S 2 )], ^l p0 = 
exp[— t 2 /(2Tp)], and G = ir/4. The plots clearly show that the idler is the time-reversed and 
phase-conjugated replica, i.e. SPC, of the signal. 




Fig. 4. Signal gain rj + 1 and idler gain rj versus G from numerical analysis compared with 
theory. See caption of Fig. for parameters used. 
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Figure |3] plots the idler gain on the logarithmic scale for a wider range of G's and two different lengths, 
obtained from the numerical analysis of the complete three- wave-mixing equations (J5J, and JBJ, with a 
single photon as the input signal, approximately emulating parametric fluorescence. For the L = 10 cm case 
the curve can be clearly separated into three regimes; for G < n/2 and moderate gain (77 ~ dB), the idler 
gain approximately follows the Fourier solution (dashed curve). For G > n/2, the system becomes unstable 
and an exponential growth (linear ramp on the logarithmic curve) is observed, until the pump is significantly 
depleted, parametric oscillation occurs and the exponential growth abruptly stops. 

150 



100 



g 50 





""0 0.5 1 1.5 2 2.5 

G 

Fig. 5. Plot of numerical idler gain 77 in dB against G for L = 10 cm (solid) and L = 1 cm 
(dash-dot), compared with the Fourier theory (dash), tan 2 (G) in dB. Three distinct regimes 
can be observed for the L = 10 cm case; the moderate gain regime where the Fourier theory 
is accurate, the unstable regime where the gain increases exponentially, and the oscillation 
regime where significant pump depletion occurs. For L = 1 cm, the medium is not long 
enough for oscillation to occur in the parameter range of interest. 

For L = 1 cm, the numerical solution departs from theory for a smaller G, and the slope of the logarithmic 
curve in the unstable regime, proportional to L, is too small to initiate oscillation in the parameter range of 
interest. 

A medium length of 10 cm may be pushing the limit of current technology. Even if one is able to fabricate 
such a long periodically-poled nonlinear crystal, the effective medium length is always limited by parasitic 
effects, such as diffraction, group- velocity dispersion, and competing third-order nonlincaritics, so it might be 
difficult to fabricate an ideal EPM device for the aforementioned purposes. For instance, in the experiment 
by Kuzucu et at, the diameter of the beam is W ~ 200 /_tm, so the characteristic diffraction length is 
<~ W 2 /Xo = 4 cm, while the characteristic group-velocity dispersion length is 20 cm according to Ref. 16, 
which are all on the order of the medium length required for mirrorless OPO. That said, techniques like 
diffusion bonding 34 can be used to increase the length of a nonlinear crystal, diffraction can be eliminated 
by waveguiding, while there exist a variety of methods to compensate for group-velocity dispersion and 
third-order nonlinearities. 35 Hence with careful engineering, fabricating an EPM device for the proposed 
applications is still a distinct possibility. 
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7. Conclusion 

In summary, it is proven that the copropagating three-wave-mixing process, with appropriate extended 
phase matching and pumped with a short second-harmonic pulse, is capable of performing spectral phase 
conjugation, parametric amplification and efficient multiphoton entanglement. The main technical challenges 
of experimental implementation seem to be the long medium length required and the control of parasitic 
effects such as diffraction, group-velocity dispersion, and competing third-order nonlincarities. However, a 
shorter proof-of-concept device has already been experimentally realized for the purposes of broadband 
second-harmonic generation 18 and coincident frequency entanglement, 17 so it is not unrealistic to expect 
that a longer device can be fabricated for the proposed applications, which should be useful for optical 
communications, signal processing, and quantum information processing. 

Theoretically, much remains to be explored. The study of parasitic effects, not considered in this paper, is 
vital for experimental realization. The analysis of the ultrashort-pump limit can be potentially generalized 
to other TWM and FWM geometries, while the quantum analysis of this limit is by no means complete. In 
conclusion, the analysis presented here should stimulate further experimental and theoretical investigations 
of a new class of parametric devices. 

The author would like to thank Prof. Demetri Psaltis for helpful discussions and a reviewer for pointing 
out Rcfs. 19-22. 
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